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ABSTRACT 

We investigated the structures on scales beyond the typical clusters of galaxies. 
These structures are crucial to understand the cosmic gravitational clustering in the 
pre-virialized stage, or quasilinear regime. Based on the multi-resolution analysis of the 
discrete wavelet transformation, we got statistical available ensembles of rc^-clusters, 
i.e. the structures on scale rd, in the range 1 < r^i < 24 Mpc for both N-body 
simulation and the IRAS 1.2 Jy galaxy survey samples. If the mass-to-light ratio 
on scales larger than clusters asymptotically reaches a constant, we found that the 
abundances and correlations of these IRAS rc^-clusters to be basically consistent with 
the predictions of the flat low-density CDM model (LCDM) and the open CDM model 
(OCDM), except the model-predicted abundance of rci=2A Mpc clusters seems to 
be higher than IRAS data. The standard CDM (SCDM) gives too much power, and 
too weak correlations on all the scales. For a given rd, the amplitude of two-point 
correlation function of r^z-clusters is increasing with their richness. However, for a 
given richness (defined by the mean separation of neighbor objects), the clustering 
strengths of both simulation and observation sample are found to be declining with 

when is larger than 3-4 Mpc. Therefore, the "universal" increase of the 
correlation amplitude with the scale of objects from galaxies, groups, to poor, rich 
clusters is broken down for structures of > 3 - 4 h^^ Mpc. Supercluster should not 
be a member of the "universal" increase family. 



Subject headings: cosmology: theory - galaxies: clusters: general - large-scale structure 
of universe 



1. Introduction 

The study of clusters of galaxies is among the most fruitful topics in cosmology. It provides 
rich information on large scale structures, significant for understanding the nature of formation and 
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evolution of cosmic structures on scales of ^ 1 - 10 Mpc. In particular, cosmological parameters 

such as the density parameter fio, the mass fluctuation amplitude as and even the cosmological 
constant 0,x can be effectively constrained by observed properties of clusters. For instance, low 
mass density cold dark matter (CDM) models, fio — — 0.3, are found to be generally favored by 
the abundance and two-point correlation function of clusters (Jing Sz Fang 1994; Eke et al. 1996; 
Viana & Liddle 1996; Carlberg et al. 1997; Bahcall, Fan & Cen 1997; Croft et al. 1997). 

However, the redshift and peculiar velocity data on scales larger than clusters, such as IRAS 
galaxy survey, yield CIq ~ 0.3 — 1.5, with a trend towards — 1 if moderate biasing is taken into 
account (e.g. Dekel 1994; Strauss & Willick 1995 for summary of the results). Moreover, structures 
on scales larger than typical clusters are pre-virialized, or depend on the quasilincar evolution of 
cosmic clustering. Namely, the larger-than-cluster (LTC) objects are in a quantitatively different 
dynamical evolution stage from the objects on scales equal to or less than Abell rich clusters 
which are virialized, or nearly virialized. Thus, one cannot estimate the clustering behavior of 
larger-than-cluster (LTC) objects by the results of clusters. Therefore, it is worth to investigate 
whether the low mass density CDM cosmogony is still favored by the abundance and two-point 
correlation function of structures on scales larger than clusters. This is the goal of this paper. 

A curious thing is that we know of some observed LTC structures, like superclusters, and we 
know the dynamics of structure formation of popular models, but at the moment have not put 
them together. This is largely because the lack of a well-defined ensemble of LTCs. Structures on 
scales larger than clusters generally are not "regular" in both configuration and velocity spaces. 
There is no weU-defined dynamical criterion to classify and identify these structures, and therefore, 
statistical approaches have not yet been carefully pursued. 

This problem can be solved by a complete decomposition. Since a complete decomposition 
does not lose information, any "irregular" objects or distributions can be precisely represented 
by the coefficients of a complete decomposition. Thus, the model-observation confrontation can 
be done with a decomposed representation if the decomposition can reasonably work for both 
observed and model-predicted distributions. Recently it has been shown that a statistical analysis 
can indeed be performed for un-virialized large scale structures by a complete multi-resolution 
analysis of the discrete wavelet transformation (DWT). It has successfully identified structures on 
scales from 1 to about 20 Mpc, and revealed the scaling behavior of LTC clustering (Xu, Fang 
& Wu 1998, hereafter Paper I). 

This paper is organized as follows. In §2 we present the abundances of IRAS rcj-clusters and 
compare with the model predictions. §3 tests the models with the IRAS two-point correlation 
functions of the LTC structures. Finally, conclusion is in §4. 



2. Abundances of dusters 
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2.1. Mass functions of rcz-dusters from simulation samples 

The model samples of mass distributions used in this paper are generated by N-body 
simulations with the P^M code developed by Y.P. Jing (Jing & Fang 1994, Jing et al. 1995). Three 
popular models of the CDM cosmogony, i.e. the standard CDM (SCDM), the flat low-density 
CDM (LCDM) models and the open CDM model (OCDM), arc employed. The Hubble constant h 
(= Hq/(WO km/scc/Mpc)), mass density Qq, cosmological constant Qx cind the normalization of 
the power spectra of the density perturbation dg are chosen to be (0.5, 1.0, 0.0, 0.62) for SCDM, 
(0.75, 0.3, 0.7, 1.0) for LCDM, and (0.75, 0.3, 0.0, 1.0) for OCDM. The simulation parameters are: 
1.) box size L = 310 Mpc; 2.) number of simulation particles Np = 64^ and 3.) the effective 
force resolution r) = 0.24 h"^ Mpc. In this case, a particle has mass of 3.14 x 10^^QMh''^MQ, 
which is small enough to resolve reliably objects with mass M > 5.5 x 10^^/i~^Mq. We produced 
8 realizations for each model. The mass of a rc/-cluster is measured within a sphere of radius rd- 

The method of finding clusters on scales rd by discrete wavelet transformation (DWT) has 
been developed in Paper I (about the technical aspects of the DWT multi-scale identification, sec 
also Fang & Thews 1998). We identified clusters on scales rd = 0.75, 1.5, 3, 6, 12, and 24/j.^^Mpc 
from the simulated samples. A result is shown in Fig. 1. It plots the top 25 massive rci-clusters for 
each scale of rd = 0.75, 1.5, 3, 6, 12, and 24 Mpc, from one realization of OCDM simulation 
with volume 310^h~^ Mpc^. 

Observationally, supercluster is defined as cluster of rich clusters (top Tab-clusters). In Fig. 1, 
some Tab-clusters are located together to form large scale clumps or filaments. These structures 
are obviously superclusters. Fig. 1 also shows, however, that some of the top 25 rd = 24 
Mpc halos contain nothing of rd < 24 Mpc massive halos. This is, although the mass and 
length scales of these rd = 24 Mpc clusters are the same as superclusters, they do not 
contain any of the top 25 clusters on smaller scales. This result indicates that the identification 
of superclusters by "cluster of rich clusters" may miss some very massive objects on scale rd = 
24 Mpc. Recently, it has been found that the dynamical properties of superclusters are not 
only determined by rich clusters, but also by numerous galaxies in the intra-cluster regions (Small 
et al. 1998). This also implies that one should directly identify LTC structures from matter or 
galaxy distributions, instead of using only rich clusters. 

With the multi-scale identification of simulation samples, we have model-predicted mass 
function of Tc^-clusters, n(> M,rd), which is the number density of rc^-clusters with mass larger 
than M. The mass functions of clusters with rd = 1.5, 3, 6 and 12 Mpc in models of the 
SCDM, LCDM and OCDM are shown in Fig. 2. The mass function of clusters on Abell scale, 
i.e. rd = rah = 1-5 Mpc, is the same as that given by friends- of -friends identification (Jing 
& Fang 1994). We have also calculated mass functions of rd = 24 Mpc clusters, which are 
not shown in Fig. 2, but will be employed below. The statistical errors of the mass functions of 
Tcj-clustcr are mainly caused by Poisson errors. The middle part of the mass function is usually 
more reliable, as all the Poisson errors are relatively small in this range. In Fig. 2, this middle 
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range corresponds roughly to 10 ^ < n(> M,rci) < 10 Mpc ^. 

Since each cluster now is characterized by two parameters M and Vd, it is inconvenient to 
define richness by mass alone. Instead, the relative richness can be described by the mean neighbor 
separation d between the clusters considered. The larger the d, the higher the richness. 

2.2. clusters from IRAS samples 

Wc have also applied Z?W^T multi- resolution analysis on real samples of the rcdshift survey of 
the Infrared Astronomical Satellite (IRAS) galaxies with flux limit of /go > 1.2 Jy (Fisher et al. 
1995). The IRAS 1.2 Jy sample contains 5313 IRAS galaxies, covering the whole sky with depth 
of about 150/i,-^ Mpc. 

The sampling rate of the IRAS sample is low. Therefore, the IRAS clusters identified on 
Abell scale Vab contain much fewer galaxies than optical clusters. This gives rise to higher Poisson 
errors. However, larger Vd clusters contain more members, and Poisson error is lower. The IRAS 
sample can give statistically significant ensemble of r^i > Tab clusters. Moreover, the IRAS sample 
consists of less early type of galaxies, and more late type. Namely, it has less galaxies in the center 
of ra5-clusters, but more field galaxies. This is rather a good point for studying structures beyond 
clusters, as late type galaxies probably are more fair tracers of mass field on larger scales. 

To minimize the radial effects, we've divided the data into redshift intervals, or the velocity 
intervals of [2500,5000], [5000,7500], [7500,10000] and [10000,12500] km s^^. We then identify 
rc/-clusters in each interval. The details of cluster-finding algorithm has been discussed in Paper I. 

It has been shown that IRAS galaxies basically trace the local large scale structures. Some 
individual clusters and superclusters from optical surveys have also been identified from the IRAS 
data (Fisher et al. 1995; Webster et al. 1997). We checked this point further by cross-identification 
between the Abell clusters and the multi-scale identified r^i clusters of the IRAS galaxies. Using 
the position data (Z, 6, and z) collected in the compilation of the Abell catalog (Abell, Corwin & 
Olowin 1989), one can find 18 rich clusters (with galaxy counts > 30) in the range [5000, 7500] km 
s-^. They are A194, A347, A400, A426, A539, A569, A779, A1367, A1656, A2870, A2877, A2911, 
A2933, A3389, A3537, A3656, A3706, A3879. 

Since the maximum positioning uncertainty of our cluster finding algorithm is about 3 
Mpc in 2-d, and 2000 km s~^ in radial, we choose the matching radius r between the projections 
of centers of the cross-identified clusters to be less than 3 Mpc, and the radial separation d of 
the centers less than 2000 km s~^. In this case, 7 out of the 18 Abell clusters are cross-identified 
with IRAS Tci = 3/i-^Mpc clusters. These are A347, A426, A539, A1367, A3389, A3656, A3706. 
The overall confidence level of this cross-identification can be estimated by a binomial distribution 
function B{k,n,p) that a random distribution would generate > k cross-identifications from n 
realizations when the chance alignment probability is known as p. When overlapping is not 
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important, p = Nirr'^d/V, where N and V are the total number of clusters and the volume of the 
redshift shell, respectively. For the above-mentioned Abell-IRAS cluster cross-identification the 
overall confidence level is Scr (~ 99.7%). 

The difference between 7 and 18 is due partially to the low sampling rate of the IRAS, i.e. the 
higher Poisson errors. It is not a surprise that the richness of individual clusters estimated from 
IRAS may have large difference from optical estimate. This cross-identification seems to show 
that the Abell clusters may be underestimated by the IRAS identification by a factor about 2.6. 
However, we also did the cross-identification in the opposite direction, i.e. to identify the 16 top 
IRAS clusters in the shell [5000, 7500] km with Abell clusters. The result is: 5 out of the 16 
IRAS clusters have their Abell counterparts. Namely, the IRAS clusters may be underestimated 
by the Abell clusters by a factor about 3. Therefore, satistically, both Abell and IRAS clusters 
trace the same density field. If we adjust the matching radius r and separation d, we would find 
that the number of Abell-IRAS cross-identified clusters increases from the 7 of 18 (3cr confidence 
level), to 13 {2a) and 15 (Icr). This shows that the Abell and IRAS clusters provide about the 
same description of the massive halos on scale rab- 

Using the 7 better cross-identified ra = 3h~^Mpc IRAS-Abell clusters, one can calculate the 
mean ratio of NiRAs/Nabell, where Njras = Ng{l.b) /<i){z)), Nabell is the number of galaxies of 
the considered Abell cluster, A'^,^(1.5) is the observed number of IRAS galaxies within radius 1.5 

Mpc around the cluster center, and (j){z) the IRAS selection function (Fisher ct al. 1995). 
We found that the median value Niras /^abeii = 7.03. On the other hand, the mass-to-number 
ratio of Abell clusters on average is Mabeii/Nabell = 0-6 x lO^^^^^M©, where M^beii is the mass of 
the Abell clusters (Bahcall and Cen 1993). Thus, the total mass of the IRAS clusters M can be 
represented as 

M = ANiRAs, (1) 
where the number-mass conversion coefficient A is 

A = 10^^-^^ Mq. (2) 

2.3. Mass function of IRAS rab-clusters 

Considering the luminosity is proportional to the number of galaxies, the number-mass 
conversion relation eq.(l) essentially is the mass-to-light ratio for IRAS Tab-clusters. Moreover, the 
mass-to-light ratio for Abell clusters is independent of the total mass of the clusters (Bahcall & 
Cen 1993). Therefore, it seems to be reasonable to employ eq.(l) to find the total mass of IRAS 
Tab-clusters. 

The total number of member galaxies, Njras, is given by Njras = ^g/^{z), where (f){z) is 
the selection function (Fisher et al. 1995). Thus, we have 



M = ANg/cj){z). 



(3) 
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Using eq.(3), the mass function of the IRAS Tab-clusters can be found from their abundance 
nnc{> Ng,rci) by 

niRAs{> M, rab) = n„c(> Ng = Mcf){z)/A, r^b). (4) 

Using the nnc{> Ng,rab) given by the IRAS galaxies in the shell [5000, 7500] km s~^, 
we calculated the mass function of the IRAS r„6-clusters nijiAs{> M,rab)- We found that 
n'lRAs{> M,rab) is in very good agreement with the mass function of optical and X-ray clusters 
within Abell radius, nop{> M,rab) (Bahcall and Cen 1992), espicially for rich clusters of M > 10^^ 

We have also calculated the mass function niRAs{> M,rab) using the IRAS data of the shells 
[2500, 5000], [7500, 10000] and [10000, 12500] km s^^ The resulted niRAs{> M,rab) is also 
basically in agreement with nop{> M,rab)- The goodness of this agreement can be desrcibed by 
the number-mass conversion coefficient A. In order to have a best fitting of niRAs{> M,rab) with 
nop{> M,rah), we find log^ = 12.07 ± 0.21 for shell [2500, 5000] km s"!, 11.93 ± 0.15 for [5000, 
7500], 11.71 ± 0.32 for [7500, 10000], 11.63 ± 0.41 for [10000, 12500]. Here the errors include both 
the Poisson errors of tijras and riop- A shows a slight decrease from lower velocity shell to higher 
velocity. This might be caused either by the uncertainty of the selection function, or by a very 
weak redshift variation of IRAS galaxies with respect to optical clusters. 

The conversion coefficients A from all the redshift shells are the same as the A of eq.(2) within 
errors, and their differences are not larger than a factor of 3. This is, in terms of mass function 
of Tab-clusters, the matter distributions traced by the IRAS galaxies is about the same as optical 
and X-ray clusters on large scales. This result seems to be a little surprising. It is known that 
the luminous IRAS galaxies with heavy dust lane do not appear at the densest regions, namely 
the clustering of infrared galaxies is weaker than that of optical galaxies at the core region of 
rich clusters. Therefore, the matter distributions represented by IRAS galaxies should be less 
clustering than that traced by optical galaxies on small scales. In other words, the mass-to-light 
ratio should actually be higher for IRAS clusters at the core region. However, the mass function 
given by eq.(4) matches well with that of the optical clusters. This is because the mass function 
eq.(4) measures statistically the mean number of collapsed massive rafo-halos over a large volume. 
It is actually a measure of clustering features on scales equal to or larger than clusters. It thus 
less depends on the physical process within the core region. In a word, the IRAS galaxies may 
cause underestimates of individual masses of some richest clusters on smaller scales, while for large 
scales, IRAS galaxies give about the same statistically features of mass field as optical objects. 



2.4. Model comparisons of abundances of IRAS rd-dusters 

Similar to eq.(4), one may derive the mass function of IRAS rc/-clusters, niRAs{> M,rci), 
from their number-count function nnc{> ^g^fci) by 

niRAs{> M, Vol) = n„c(> Ng = M(t>{z)/A, ra). ( 
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The problem is yet how to determine A at scale rd > Tab- 

Prom eq.(l), A essentially is the mass-to-light ratio of the system on scale Vah- The derived 
value of A for rab-clusters corresponds to a median mass-to-light ratio of M/L ~ SOOiblOO h{M/L)Q. 
As mentioned above, rat-clusters are dominated by elliptical galaxies, while the larger structures 
have higher fraction of spiral galaxies. The mass-to-light ratio or A on scales rd > Tab may be 
smaller than Tab- The observed Virgo cluster infall motion, however, indicates that M/L ratio 
doesn't show a decrease with scale, while appears to flatten and remain approximately constant 
on scales larger than 0.2 h"^ Mpc (Bahcall, Lubin & Dorman 1995). This result is strengthened 
by a direct measure of M/L of the Corona Borealis supercluster (Small et al. 1998). It gives 
M/L ~ 560 h{M/L)Q. Namely, the M/L ratio on scale ~ 20 Mpc is not lower, but even 
higher than that of rich clusters by a factor about 1.9. This implies that the M/L ratio is weakly 
scale-dependent from scales of 1 to about 20 Mpc. The weak dependence of M/L on scale is 
also supported by the following fact: the segregation between elliptical and spiral galaxies of APM 
Bright Galaxies Catalogue (BGC), as well as the ratio between the second order bias factors of 
these APM-BGC elliptical and spiral galaxies, is found to be basically scale-independent from 1 
to 20 Mpc (Loveday et al. 1995, Fang, Deng & Xia, 1998). It would be reasonable to assume 
that the elliptical-spiral mix doesn't cause a serious scale-dependence of the mass-to-light ratio on 
large scales, and A is scale-independent within a factor of 3. 

Using the A for each shell (§2.3), we calculated the mass functions of IRAS clusters with 
ra = 1.5,3,6,12 and 24h~^ Mpc for data of the velocity shells [2500, 5000], [5000, 7500], [7500, 
10000] and [10000, 12500] km s~^. The errors of the IRAS cluster mass functions are estimated 
by the same way as that of optical clusters. The vertical error bars in Figs. 3-6 are the Poisson 
errors of counting the clusters in the limited volume of a shell. The horizontal error bars are the 
Poisson errors of mass estimation of each clusters. This Poisson error is generally smaller for 
larger rd, as more galaxies are contained in larger rd clusters. Therefore, at larger rd the IRAS 
data can set effective constraints. 

Since the ratio between mass functions at rd and rab is independent of the value A, Figs. 3 
- 6 essentially are tests of the rcj-dependence of mass functions of rd clusters. Despite the errors 
of these figures are large. Figs. 3 and 4 have shown that all n(> M,rd) in the range of rd < 12 
h"^ Mpc and in the two shells [2500, 5000] and [5000, 7500] km are in good agreement with 
the predictions of the models of LCDM and OCDM. The results of 24 h~^ Mpc in these two shells 
(which have not directly been shown in Figs 3 and 4) are also basically consistent with LCDM and 
OCDM. The dynamical status of ra5-clusters is very different from rd >> rab clusters. The former 
are most likely virialized, the later are not. Therefore, Figs. 3 and 4 indicate that the models of 
OCDM and LCDM are still in good shape for clustering in pre-virialization regimes. 

For shells of [7500, 10000] and [10000, 12500] km s-^ Figs. 5 and 6 show that the LCDM 
and OCDM can fit with the data of rd < 6 Mpc. Yet, the IRAS mass functions on scales 
of rd = 12 and 24 Mpc are lower than the predictions of LCDM and OCDM. The deficit of 
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IRAS rd = 24 Mpc clusters may result from the poorness of the current data at > 7500 km 

s~^. Since the overdensity of the large scale objects is relatively low, their detection will be poorer 
than small scales. Moreover, not all predicted r^i = 12, 24 Mpc clusters contain top rich 
Tab-clusters (Fig.l), and therefore, some of these structures will be overlooked if the completeness 
of samples is not high enough. 

Figs. 3-6 clearly show a systematic and significant lower of the IRAS mass functions in all 
velocity shells than the predictions of SCDM. SCDM gives too many massive clusters on all scales 
of rd < 24 h"^ Mpc. It is interesting to note that if all IRAS mass functions are shifted along the 
log M axis by a factor of about 0.60, the SCDM model will be able to survive. The shifting factor 
0.60 requires A to be larger than our derived value by a factor lO*^'^ ~ 4, and the corresponding 
ratio M/L should be as high as 4 x 300 ~ 1200h{M/L)Q. Thus, the SCDM can not be saved if 
the uncertainty of A is less than a factor of 4. 

3. Two Point Correlation Functions of rcrClusters 

3.1. d- dependence of correlations 

It has been known for a long time that the two-point cluster correlation function ^cc{i~) follows 
the same power law as galaxies, i.e. ^cc('') = i'>'/fo)^''' and 7 ~ 1.8, but has much larger correlation 
length To than galaxies. Moreover, the correlation length rg defined by ^(rg) = 1 is found to be 
systematically increasing with the cluster richness. An empirically "universal" scaling vq = OAd, 
d being the mean neighbor separation as employed to describe the cluster richness (§2.1), was 
proposed to fit with the correlations of galaxies, QSOs, group of galaxies, X-ray and Abell clusters 
and even superclusters (e.g. Bahcall & Soncira 1983; Bahcall &; West 1992 and references therein). 
Because the sizes of these objects, from galaxies to superclusters, are systematically increasing by 
at least two orders, it is unclear whether various objects with very different scales can really be 
fitted with one "universal" scaling. 

This problem can be solved by studying the two-point correlations of clusters with different 
rd- Using samples of rc^-dusters from N-body simulation, we found that the two-point correlation 
functions of rcrclusters can be well fitted with a power law ^(r) = {r/ro)~'^ for all considered rd 
and d. The index 7 is in the range of 1.8 - 2.0 for the LCDM and OCDM, and 1.9-2.1 for the 
SCDM. No indication shows that the index 7 would significantly vary with either rd or d. 

The ro-d relations for rcrdusters are calculated. The results of the OCDM are plotted in Fig. 
7. This figure shows that the ro-d lines for rd = 0.75, 1.5 and 3 Mpc almost coincide. This 
is, in the range of Vd < 3 Mpc, the r^-d relation is "universal", independent of the size of 
structures. The model SCDM and LCDM show similar "universal" ro-d relations. 

However, the r^-d relation is substantially dependent on rd if rd > 3 Mpc. Despite 
the correlation length ro is still increasing with richness d, the lines ro — d are lower for larger 
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rd- Therefore, superclusters are not a member of the universal rQ-d relation of < 3 /i Mpc 
clusters. 

For clusters of Vd larger than 12 hT^ Mpc, their ro is comparable to and even less than rd- 

In this case, ro is actually determined by an extrapolation of the two-point correlation function 
^(r) from range r > Vd to r < rd, and solving ^(^o) = 1. This method might be susceptible to 
large error. Nevertheless, it is reliable that lines of ro — d are significantly lower than the ro — ci of 
smaller rd clusters. 

In redshift space the correlation length will increase by no more than 2h^^ Mpc, and the 
index 7 becomes slightly larger. Our overall conclusions on correlation functions will not change 
due to these slight differences. 

3.2. rcrdependence of correlations 

Fig.8 plots the correlation length ro as a function of Vd for a given d. This figure, actually, 
is an alternative representation of Fig.7. It shows that for a given d the correlation length ro is 
almost invariant for rd equal to or less than 3 Mpc, but slowly decreasing with Vd when rd is 
larger than about 4 Mpc. This is a common feature for all the three models. 

Theoretically, the decrease of rg with rd is expected. Fig. 1 shows that most halos of 
rd = 0.75, 1.5 or 3 /i^^ Mpc coincide to form a 0.75-1.5-3 clouds-in-clouds. This is because most 
halos of rd < 2 h^^ Mpc are virialized. Therefore, the detection of rc/-clusters with rd < 3 
Mpc will find almost the same objects. In other words, the ensembles of rd = 0.75, 1.5 or 3 
Mpc clusters actually consist of the same objects, and therefore, they have the same correlations, 
and the universal ro-d relations. However, on scales larger than 4 Mpc, more and more halos 
are less developed, and their correlations are lower. Therefore, the "universal" increase of the 
correlation amplitudes with scale from galaxies, groups to clusters should be broken down for rd > 
3-4 Mpc clusters. Superclusters are not members of the "universal" scaling family. 

3.3. Correlation lengths of IRAS clusters 

With the samples of the IRAS rcj-clusters, wc calculated the correlation length rg by fitting 
the observed data with ^(r) = (r/ro)~^'^. The correlation lengths ro are measured from data of 
different velocity shells. The results of ro vs. d are plotted in Fig. 9. The data points shown in 
Fig. 9 are binned in d. The errors are given by the variance of the correlation strengths at each 
bin of d. For rd = fab, the values of ro found from the IRAS samples arc in agreement with those 
from APM sample (Bahcall & Cen 1992; Croft et al. 1997). Fig. 9 shows, obviously, that the r^-d 
relation from IRAS samples of rd = 1-5 to 12 Mpc favors the models of OCDM and LCDM. 
The SCDM-predicted rQ-d curves are systematically lower than the IRAS results for all considered 
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d. 

We also tested the models by the rcr dependence of Vq. Fig. 10 gives the ro-Vd relations from 
the simulation samples and IRAS data. The errors of the IRAS data are also large. But the 

observed data show once again that the LCDM and OCDM are basically consistent with the IRAS 
data, while the SCDM-prcdictcd tq is too low to match with the observation. More importantly, 
it is clear to see the slow decrease of vq with the increase of rd in all models. The ro decrease is 
an indicator of the transform from quasilinear to non-linear regimes. 

4. Conclusions and discussions 

Structures on scales equal to and larger than typical clusters have been identified by the 
Z)M^r decomposition. This method is shown to be able to systematically and uniformly identify 
structures on various scales, especially those of larger than the size of typical clusters. Taking 
this advantage, one can produce statistically available ensembles of rcj-dusters from simulation 
samples and observational data. 

To describe the dynamics of the cosmic clustering on scales beyond clusters, we introduced 
several new statistical measures: the bivariant abundance function n{> M,rci), the re^-dcpcndence 
of abundance, the rcrdependence of correlations, etc. These statistics are useful to reveal the 
features of pre-virialized regime of cosmic gravitational clustering, and to discriminate among 
models of structure formation, especially on large scales. 

The identification of LTCs is effective to extract information from sparsely distributed data 
like the IRAS galaxies. The IRAS 1.2 Jy galaxies have low sampling rate compared to optical 
galaxies and are not found at cluster centers. They are not very effective to be employed for 
identification of traditional raj-clusters. However, the IRAS clusters with Vd > 3 h^^ Mpc can 
already set effective constraints on models, if the mass-to-light ratio for structures traced by 
IRAS galaxies is assumed to asymptotically reach a constant on scales larger than Abcll radius. 
Therefore, IRAS galaxies can provide the same information of large scale structures as optical and 
X-ray galaxies, if the considered scales are larger than the scale of morphological segregation of 
IRAS galaxies that no dusty spirals like IRAS sources exist at the immediate neigbhourhood of 
densest regions. 

The multi-resolution analysis of the IRAS sample concludes that the abundances n{> M,rci) 
and correlations of objects on scales rd = 1 to 24 Mpc can be described by LCDM and 
OCDM, while the SCDM is in difficulty. Namely, the abundances and two-point correlation 
functions of structures on scales as large as about 20 Mpc favor a low-density universe if M/L 
ratio is assumed to be flat over these scales. 
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Fig. 1. — A projected distribution of rc;-clusters identified from a realization of OGDM N-body 
simulation within box 310 Mpc^. It shows the top 25 of the richest clusters for each ra- The 
legends of the symbols are: plus for 0.75, diamond for 1.5, circle for 3, triangle for 6, square for 
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12 and giant circle for 24 h ^ Mpc. The linear sizes of symbols are roughly equal to their rd, but 
slightly amplified for small rd for clarity. 

Fig. 2. — Mass functions of clusters identified with radii Vd = 1.5, 3.0, 6.0 and 12 Mpc 
for SCDM, LCDM and OCDM models at redshift z = 0. n(> M,rd) is the number density of 
Tcrclusters with mass larger than M. M is in unit of h"^ Mq. 

Fig. 3. — Mass functions of IRAS rcrdusters. n is in unit of Mpc~^. The horizontal errors are 
Poisson errors from the number counting of galaxies in a given cluster, and the vertical errors are 
the Poisson errors from counting the clusters. The model predictions of SCDM, LCDM, OCDM 
are shown as dashed lines, thin solid lines, and thick solid lines, respectively. Observed data are 
from the IRAS sample in the velocity shell [2500, 5000] km s~^. 

Fig. 4. — The same as Fig. 3, but the data are from shell [5000, 7500] km s~^. 

Fig. 5. — The same as Fig. 3, but the data are from shell [7500, 10000] km s~^. 

Fig. 6. — The same as Fig. 3, but the data are from shell [10000, 12500] km s~^. 

Fig. 7. — The rQ-d relations for clusters with different r^i in the OCDM model. The error bars of 
the curves of rd = 6 and 12 Mpc are from the scattering among realizations of simulation. 

Fig. 8.— The ro-rd relations at different richness for models of OCDM, LCDM, SCDM. 

Fig. 9. — The ro-d relations for rd = 1.5, 3, 6, I2h~^ Mpc. The dashed, thin solid, and thick solid 
lines are, respectively, the simulation results of models SCDM, LCDM and OCDM. The observed 
data are from the averages over the velocity shells of the IRAS galaxies. 

Fig. 10. — The ro-rd relations for d = 25, 50, 70 h"^ Mpc. The dashed, thin solid, and thick solid 
lines are, respectively, the simulation results of models SCDM, LCDM and OCDM. The observed 
data are from the averages over the velocity shells of the IRAS galaxies. 
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